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ABSTRACT 


A  method  to  infer  sediment  concentrations  from  acoustic  backscatter  levels  has  been 
implemented  for  a  prototype  acoustic  sediment  flux  probe.  This  requited  calibrations  of  the 
acoustic  transceiver  systems  and  direct  measurements  of  the  system  response  to  typical 
sediment  sue  distributions  over  a  wide  range  of  concentrations. 

The  sensitivities  of  the  1.3  MHz  and  5.2  MHz  transducers  of  an  acoustic  sediment  probe 
were  measured  using  the  backscatter  amplitudes  from  stainless  steel  wires  of  four  different 
radii.  A  two-frequency  inversion  algorithm  estimating  the  geometric  mean  radius  and 
variance  about  the  mean  of  suspended  sediment  with  an  assumed  lognormal  size  distribution 
and  the  sediment  mass  concentration  was  developed  and  tested  in  a  laboratory  setting. 

The  sensitivity  ofthe  5.2  MHz  transducer  for  each  wire  was  consistent  within  10%.  The 
sensitivity  of  the  1.3  MHz  transducer  calculated  for  each  wire  was  not  consistent  and  is 
believed  to  be  due  to  a  strong  angular  dependence  between  die  wire  orientation  and  foe  1 .3 
MHz  transducer  face.  The  sensitivities  of  each  of  the  four  transducers  were,  however, 
inferred  from  the  system’s  response  to  the  controlled  sediment  backscatter  measurements. 
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L  INTRODUCTION 


The  coastal  region  of  the  United  States  is  of  high  economic,  recreational  and  strategic 
value.  Ports  and  harbors  are  important  for  the  import  of  foreign  goods  and  the  export  of 
agricultural  and  industrial  products.  The  coastal  region  is  subject  to  many  natural  forces, 
often  resulting  in  disastrous  consequences.  Although  storms  have  important  effects  on 
coastal  property,  the  continual  exposure  of  the  nearshore  zone  to  waves  and  currents  can 
have  cumulative  effects  which  can  be  significant  over  a  long  period  of  time.  Depending  on 
die  wave  dynamics  and  nearshore  topography,  alongshore  littoral  drift  can  cause  erosion, 
transport,  and  deposition  of  vast  amounts  of  sediment  The  potential  for  extensive  damage 
to  shoreline  homes,  real  estate,  and  harbor  structures  exist  in  both  die  short  and  long  term. 
Prediction  of  sediment  transport  rates  require  that  die  amount  of  suspended  sediment  flux 
be  known  undo'  a  variety  of  wave  conditions.  As  deep  water  ocean  waves  and  locally 
generated  wind  waves  propagate  toward  the  coast  their  energy  is  transformed  to  higher  and 
tower  frequencies  in  nonlinear  interactions  under  die  influence  of  die  underlying  topography. 
The  nodnear  fluid  motions  cause  transport  of  die  sediment  at  die  seabed,  potentially  leading 
to  significant  changes  in  die  local  topography.  Much  of  die  incident  wave  energy  is 
dissipated  by  small-scale  turbulent  motions  which  strongly  couple  larger-scalc  fluid  motion 
through  die  water  column  to  die  sediment  bed.  Small-scale  turbulent  kinetic  energy  is 
produced  at  both  the  surface  by  non-linear  wave  interactions  and  wave  breaking,  and  at  die 
bottom  boundary  due  to  shear  stress.  While  the  bottom  boundary  layer  is  typically  only  a 
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few  centimeters  thick  under  the  oscillatory  flows  found  in  the  nearshore  region,  under  more 
energetic  wave  breaking  conditions  the  turbulent  boundary  layer  can  occupy  the  entire  water 
column.  Consequently  die  strong,  episodic,  multi-scaled  forcing  in  the  nearshore  region 
presents  unique  problems  in  our  understanding  sediment  transport  processes. 

Models  and  measurement  of  sediment  transport  in  the  nearshore  must  take  into  account 
not  only  die  strong  irrotational  wave  motion  above  the  sediment  bed,  but  also  the  much 
weaker  turbulent  motions  generated  by  the  fluid  flow.  Development  and  verification  of 
models  with  in-situ  data  requires  high  spatial  and  temporal  sampling  throughout  die  water 
column.  Ideally,  sediment  fluxes  should  be  measured  to  the  smallest  turbulent  scales 
throughout  the  water  column  to  investigate  the  episodic,  non-linear  physical  processes  which 
produce  net  sediment  fluxes  in  the  nearshore. 

Methods  of  sediment  concentration  and  size  distribution  measurements  range  from  die 
physical  collection  of  a  volume  of  suspended  sediment  and  the  medium  of  suspension,  to 
non-invasivc  methods  baaed  on  the  scatter  of  incident  acoustic  or  light  energy  by  suspended 
scatterers.  Siphoning  of  a  sample  using  suction  involves  the  removal  of  a  volume  of 
sediment  and  medium  from  die  flow  without  erroneously  sampling  the  sediment  bed.  This 
method  is  restricted  by  its  sensitivity  to  background  currents,  particle  size,  air  entrainment, 
and  die  direction  of  suction  in  relation  to  the  direction  of  flow  (Bosnian  et  al,  1987). 
Samples  are  typically  siphoned  in  a  direction  normal  to  the  flow  with  a  suction  velocity 
which  exceeds  the  ambient  flow  velocity  by  more  than  three  times  (Bosnian  et  al,  1987),  but 
samples  can  also  be  taken  along  the  axis  of  die  flow  with  a  suction  velocity  which  exceeds 
die  maximum  settling  velocity  of  die  particles  under  investigation  (Hay,  1991).  WhOe  these 
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methods  provide  a  direct  sampling  of  the  concentration,  they  clearly  perturb  the  flow  and  are 
limited  in  their  volume  resolution,  and  extremely  limited  in  temporal  resolution. 

Transmissometers  determine  die  extinction  characteristics  of  the  medium  by  measuring 
the  attenuation  of  light  along  a  short  light  path  through  the  medium.  They  provide  a  well 
controlled  sample  volume  but  have  limited  dynamic  range.  Empirical  models  are  required 
to  relate  die  beam  attenuation  to  suspended  mass. 

Optical  backscatter  instruments  measure  turbidity  and  suspended  particle  concentrations 
by  detecting  scattered  infrared  radiation.  The  response  of  an  optical  backscatter  sensor 
depends  on  the  sue,  composition,  and  shape  of  suspended  particles.  This  requires 
calibration  with  suspensions  from  the  waters  to  be  investigated  (D  &  A  Instrument 
Company,  1991).  OBS  instruments  are  wefl  suited  to  the  measurement  of  suspended 
concentrations  ith  their  wide  dynamic  range,  linear  response  to  increasing  concentrations, 
and  insensitivity  to  bubbles  and  organic  matter.  Their  primary  limitation  for  small-scale 
studies  is  that  the  sample  volume  depends  on  the  concentration  of  the  suspended  sediments 
being  measured.  This  leads  to  a  poorly  defined  sample  volume,  and  so  are  mote  suited  to 
die  measurement  of  bulk  sediment  concentrations. 

The  use  of  acoustic  backscatter  instruments  is  a  promising  area  of  research.  Small  and 
robustly  engineered  acoustic  backscatter  instruments  are  well  suited  to  sediment  flux 
measurements.  Unlike  siphon  sampling,  die  acoustic  instrument  can  be  deployed  well  above 
die  boundary  hycr  (typically  1  meter),  directed  downward  to  measure  the  turbulent  boundary 
layer  properties  with  non-mtxusive  acoustic  beams.  Current  implementations  of  this 
technology  have  high  spatial  [0(1  cm)]  and  temporal  resolution  [0(0.1  a)]  providing  a 
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method  for  sampling  the  small-scale  concentration  levels  in  the  boundary  layer  (Hay,  1983). 
Sampling  (he  sediment  along  the  length  of  die  acoustic  beam  in  range  bins  is  accomplished 
by  rapidly  sampling  received  acoustic  energy  levels  from  short  duration  transmitted  acoustic 
pulses.  Inversion  techniques  are  used  to  convert  the  acoustic  backscatter  levels  to  particle 
concentration  and  size  distribution  parameters.  This  allows  measurements  of  sediment 
concentration  properties  spanning  the  depth  of  the  boundary  layer.  The  inversion  technique 
developed  in  this  study  is  based  on  previous  work  (Hay  and  Sheng,  1992)  utilizing  multiple 
acoustic  frequencies. 

The  acoustic  backscatter  instrument  utilized  in  this  study  is  die  Coherent  Acoustic 
Sediment  Probe  (CASP)  developed  under  the  direction  of  T.  P.  Stanton  at  die  Department 
of  Oceanography,  Naval  Postgraduate  School  The  CASP  is  a  two-frequency,  dual  mode 
underwater  instrument  package  originally  designed  for  the  study  of  surface  boundary  layer 
turbulence.  It  uses  three  5.292  MHz  transducers  and  a  single  1.323  MHz  transducer  to 
determine  three  components  of  velocity,  suspended  sediment  concentration  profiles,  and 
sediment  size  distribution  using  a  multi-frequency  approach  described  in  die  following 
chapter.  By  combining  die  three  velocity  components  measured  by  the  CASP  and  the 
sodhnent  concentration  estimated  from  an  inversion  of  the  volume  backscattered  strengths, 
three  components  of  the  sediment  flux  can  be  determined  down  to  centimeter  scales. 

In  this  study,  die  absolute  sensitivity  of  the  CASP  acoustic  backscatter  system  was 
determined  using  a  fixed  wire  backscatter  calibration  technique  following  die  theoretical 
development  by  Sheng  and  Hay  (1993).  A  verification  of  the  two  frequency  inversion 
algorithm  was  made  with  the  CASP  using  controlled  concentrations  of  sediments  collected 


4 


from  local  beaches.  Acoustic  backacalter  data  were  collected  for  sub-octavc  concentration! 
up  to  2. 1  g/htcr  using  sorted  sediment  samples  of  known  lognormal  size  distributions.  A 
homogeneous  sediment  suspension  of  pre-determined  mass  concentrations  was  achieved 
using  a  uniquely  designed  test  vessel.  The  test  vessel  provided  a  well  mixed  test  volume 
with  an  acoustically  transparent  window  and  no  interior  level  surfaces  which  could  allow 
sediment  to  settle  out  during  die  trials.  The  design  also  avoided  the  generation  of  turbulent 
vortices  and  bubble  entrainment  from  the  air-water  boundary. 

The  background  theory  toed  to  determine  the  sensitivity  of  die  CASP  acoustic 
transmittor/ieceiver  and  the  sediment  inversion  technique  are  presented  in  Chapter  II.  The 
experimental  methods  used  to  determine  die  CASP  sensitivity  and  test  die  inversion 
technique  are  described  in  Chapter  HI.  Results  are  discussed  in  Chapter  IV,  and  are 
followed  by  conclusions  and  recommendations  in  Chapter  V. 
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IL  BACKGROUND  THEORY 


A.  COHERENT  ACOUSTIC  SEDIMENT  PROBE  (CASP) 

The  CASP  measures  high  temporal  and  spatial  resolution  acoustic  backscatter  information 
to  determine  the  suspended  mass,  M,  and  velocity,  V,  which  are  used  to  calculate  the  total 
sediment  flux  vector,  F~M  •  V(xy?,t) .  A  theoretical  inversion  technique  is  required  to 
convert  the  hacfcacatte  acoustic  pressure  measurements  from  the  CASP  into  sediment  size 
distributions  and  mass  concentretions.  The  amplitude  of  backscattered  acoustic  energy  at 
a  given  frequency  depends  on  the  number  of  scatterers,  their  size,  and  their  acoustic 
impedance  in  die  sample  volume.  For  simple  particle  size  distributions,  the  size  versus 
frequency  dependence  of  the  scattering  cross  sections  can  be  used  to  estimate  both  the 
scatterer  size  and  mass  concentration  by  taking  backscatter  measurements  at  multiple 
frequencies.  The  multi-frequency  approach  requires  knowledge  only  of  the  overall 
sensitivities  of  the  acoustic  instruments  and  of  the  frequency/size  dependence  of  the 
scattering  cross  section  of  the  particles. 

An  underwater  acoustic  backscatter  instrument  requires  a  transducer  to  convert  transmitter 
electrical  energy  into  acoustic  energy,  and  to  convert  received  acoustic  energy  into  electrical 
energy.  In  the  CASP,  a  transducer  of  known  dimensions  is  electrically  excited  to  produce 
a  short,  high  frequency  acoustic  pulse  which  propagates  through  the  water.  Backscattered 
energy  may  be  received  at  die  same  transducer  (monostatic  mode),  or  at  a  transducer 
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separated  in  space  from  the  transmitter  (bistatic  mode).  This  theoretical  treatment  focuses 
on  the  monostotic  mode-the  bistatic  mode  is  an  extension  using  a  different  form  factor. 

Consider  a  pulsed,  narrow-beam  transducer  shown  in  Figure  2. 1 .  The  detected  volume 
is  in  the  tor  field,  defined  by  r>  *a*/X  (Clay  and  Medwin,  1977),  where  r  is  the  distance 
from  the  transducer  to  the  detected  volume,  a0  is  the  transducer  radius,  %  is  the  attenuation 
coefficient  in  water,  a,  is  the  attenuation  due  to  the  suspended  particles,  and  A.  is  the 
acoustic  wavelength.  The  sample  volume  of  each  acoustic  beam  is  defined  across-beam  by 
the  -3  dB  beam  width,  2  p0,  and  radially  by  the  width  of  the  transmitted  pulse  length,  cx/2, 
where  c  is  the  sound  speed  in  water,  and  t  is  the  duration  of  the  transmitted  pulse.  The 
detected  volume  is  the  volumetric  intersection  of  die  cylindrical  sample  volumes  of  die 
transmitter  and  receiver  acoustic  beams.  The  transducer  detects  backscattered  acoustic 
pressure,  /^.returned  from  particles  in  the  transmitted  beam  path.  The  particles  are  assumed 
to  be  randomly  and  homogeneously  distributed  across  the  detected  volume  and  the 
concentration  is  assumed  to  be  low  enough  that  multiple  scattering  can  be  ignored.  This 
assumption  is  valid  when  die  mass  concentration  is  less  than  30  g/L  for  quartz  particles 
(Varadan  et  al,  1983).  The  maximum  mass  concentration  in  the  series  of  calibration  runs 
conducted  for  this  study  was  2.08  gfl  (32  g/15.4  liter),  although  Held  measurements  could 
well  reach  strong  scattering  values. 
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Figure  2. 1  Geometry  of  a  pulsed  monostatic  acoustic  system  and  generalized  directivity 
pattern  showing  die  main  and  side  lobes  of  die  acoustic  beam.  T  represents  the  transducer. 
The  remaining  symbols  are  defined  in  the  text. 
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ft.  TRANSDUCER  DIRECTIVITY 


The  ftrfidd  directivity,  D,  for  a  circular  piston  transducer  of  radius  %  is  given  by 

2Jt(*cfl0«nP) 

°  kca  0«inp  (2.1) 

(Clay  and  Medwin  1977)  where  J,  is  the  cylindrical  Bessel  function  of  order  1,  and 
^  (  =2*  A  ),  the  compression  wave  number.  Table  1  contains  characteristics  of  die 
GASP  transducers  for  a  sound  speed  c  *  1480  m/s  and  a  pulse  duration  t  —  2.26  x  10*3  sec. 
Tire  theoretical  directivity  pattens  for  the  1.323  MHz  and  S.292  MHz  acoustic  beams  are 
shown  in  Figure  2.2,  in  terms  of  the  spreading  angle  p. 


TABLE  1  CASP  MONOSTATIC  ACOUSTIC  BEAM  CHARACTERISTICS 


Freq. 

Transducer 

Wave- 

Faifield 

Half  Power 

Pulse 

Radius 

length 

Range 

Beam  width 

Length 

(MHz) 

(cm) 

(cm) 

(cm) 

(deg) 

(cm) 

(cm) 

1.323 

1.27 

0.112 

45.2 

2.62 

1.14 

1.68 

5.292 

0.317 

0.028 

11.3 

2.62 

1.14 

1.68 

C  BACKSCATTER  RESPONSE  TO  LOGNORMAL  SIZE  DISTRIBUTION 
The  pressure  amplitude  of  incident  spherical  sound  waves  in  tire  detected  volume  can  be 
expressed  as  (Hay  and  Sheng,  1992) 
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20lo«D)(db)  20  109(D)  (db) 


Directivity:  1.323  MHz 


Figure  2.2  CASP 1 J23  MHz  and  5.292  MHz  transducer  directivity  patterns  with  angles  of 
the  first  null  point  and  the  half  power  bandwidth  for  each  frequency  shown. 
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Pi 


P*r+  exp 


h 


-a9r  -  I  a, dr 


(2.2) 


where  />„  is  the  on  axis  pressure  amplitude  at  rm,  the  reference  distance;  r  is  the  distance 
from  the  transducer  to  the  detected  volume  with  r«r0,  is  the  attenuation  coefficient  in 
dean  under,  is  the  attenuation  due  to  die  suspended  particles.  The  overhar  denotes  the 
average  over  the  size  distribution.  D  is  the  transducer  directivity. 

Representing  the  particles  as  solid  spheres  of  equivalent  size  and  letting  The  the  output 
voltage  of  the  receiver,  it  can  be  shown  that  (Hay  1991) 


»(<*),;£■(*) 
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■wh(C) 


exp 


(2.3) 


where 


C  sct(«#+a, ) 
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2n(a)da 
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a  3  n(a)  da 
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where  S  is  the  overall  system  sensitivity  constant,  1/r  is  the  radial  spreading  term,  a  is  the 
equivalent  raefcus  of  the  particle,  n(a)  is  the  size  spectral  density,  M  is  the  mass  concentration 
of  scatterers,  and  ^  is  the  grain  density  of  the  particle.  {  is  the  combined  attenuation  of 
acoustic  energy  due  to  die  attenuation  effects  of  the  medium  and  the  scanners  in  the  acoustic 
beam.  ^  [<<*),.£(*)]  is  the  response  function  due  to  the  size  distribution  and  far-field 
backscatter  form  factor,  fjj j$  ,  of  the  particle.  Here,  x  is  given  by  x  *  k*a,  where  k.  is  the 
compression  wave  number.  The  term  |«nh(C)/C  j1*  corrects  for  attenuation  across  tbe 
detected  volume  (Hay  and  Sheng  1992). 

Most  sediment  distributions  in  the  nearshore  can  be  represented  by  a  lognormal 
distribution  of  the  form 


*(a)da  =  — — - exp 

^*Inof 


(too  -limy  f 
2(lnof)s 


(2.6) 


whiA  is  characterired  by  the  parameters  the  geometric  mean  radius  of  tbe  particles,  and 


(toof)a,  the  variance  of  In  a  about  its  mean.  For  example.  Figure  2.3  shows  a  best  fit 
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comparison  between  Hit  sue  spectral  density  for  sedimcutoolleclBdat  Monterey  Wharf  and 
a  lognormal  distribution  with  a  mean  geometric  radius  of  82  microns  and  og  -  1.22.  The 
symbob  (*)  represent  the  number  of  grains  in  each  sieve  size  range  assuming  each  sediment 
grain  has  a  diameter  equal  to  the  mean  size  of  each  sieve  range.  The  distributions  for 
Aakanar  State  Beach  and  Naval  Postgraduate  School  Beach  are  similar  with  values  of  og 
-  1.22,  -  135  microns  and  <jg  -  1.22,  =  94  microns  respectively.  This  single  mode, 

lognormal  distribution  is  a  widely  observed  sediment  size  distribution.  However,  multi¬ 
modal  distributions  can  occur,  requiring  a  more  elaborate  parameterization  using  a  larger 
wanber  of  frequencies  and  a  diffcro.  form  factor  than  the  two  frequency  approach  described 
here.  The  form  factor  for  uniformly-sized  particles,  assuming  the  particles  have  a  lognormal 
size  distribution,  is  plotted  in  Figure  2.4  as  determined  empirically  by  Hay  (1991)  and  is 
expressed 


* 


1  +1.25X4 
1+x4 


(2.7) 


where  |  is  die  theoretical  form  factor  for  a  rigid,  mobile,  spherical  particle  with 

the  density  of  quartz  (Hay  and  Mercer,  1985)  and  x  =  ,  where  <|  is  the  mean  geometric 
radius  of  the  particles.  The  dependence  of  the  form  factor  on  the  frequency  of  the  acoustic 


energy  and  die  size  of  the  particles  allows  the  use  of  two  different  frequencies  to  estimate 
the  size  of  the  particles.  The  theoretical  values  of  die  response  function. 


fln(fl),  l/LOOl  ],  for  natural  sand  particles  can  be  computed  assuming  a  lognormal  particle 
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I  (X)j{ 


Figure  2.4  Empirically  determined  form  factor  for  uniformly-sized,  spherical  particles 
assuming  a  lognormal  particle  suse  distribution. 
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size  distribution  and  the  form  fetor  for  a  lognormal  size  distribution,  [fj.x)\.  Figure  2.5 
shows  the  form  of  the  response  fraction,  given  by  Equation  (2.8)  below,  for  each  frequency 
and  values  of  of  ranging  from  1.1  to  1.5. 

After  substituting  the  tognormal  distribution  into  Equation  (2.S),  F[  n(a),  l/.OOi )  is  now 
dependent  only  on  a*,  op  and  1^.  By  letting 


-1/  2 


(lnof  f  jdlna 


J a3  exp|~(ln  a  -  In  ag )J/2  (In  of )*  j<fln  a 


i  a 


(2.8) 


Assuming  die  scattering  attenuation  is  small,  or  has  been  completely  corrected  for,  die 
backscatter  for  die  two  frequencies  can  be  expressed  by 

Kc=Si(M/ps?”FVCi,og)  i-  1,5 

where  the  subscript  i  represents  die  1.3  and  5.2  MHz  frequencies.  The  ratio  of  the 
relationship  in  Equation  (2.9)  for  each  frequency  can  be  expressed  as  (Hay  and  Shcng  1992) 

F(XS ,  Og)  =  ^<3  (2.10) 

F(Xitog)  S,  KC] 

which  eliminates  the  dependence  on  M,the  mass  concentration.  The  right-hand  side  can  be 
evaluated  from  measurements  of  the  sensitivity  S,  and  the  digital  output,  a;  ,  from  a 
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Figure  2.S  Response  functions  for  1.323  MHz  and  5.292  MHz,  assuming  a  lognormal 
particle  size  distribution,  as  a  function  of  the  particle  mean  geometric  radius  for  values  of 
of  =  1.1, 1.2, 1.3, 1.4, 1.5. 
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lognormal  sediment  sample.  The  left-hand  side  of  the  equation  can  be  evaluated  using 
Rqurton  (2.7)  using  the  expression  for  i£(x)|.  Figure  2.6  shows  the  theoretical  form  of  the 
mho  FIX, .  og)/F(Xt ,  ag).  Using  only  two  frequencies,  it  is  not  possible  to  determine  both 
ftg  and  op  so  the  value  of  the  lognormal  distribution  variance  parameter  og  must  be 
prescribed.  Theoretically  the  value  of  o  gcan  be  determined  in  methods  using  more  than  two 
frequencies,  bid  the  inherent  statistical  noise  of  the  backscatter  measurements  is  large 
enough  to  make  the  og  estimates  vary  widely,  so  og  must  be  prescribed  (Hay  and  Sheag, 
1992).  Based  on  the  values  of  o  g  determined  for  the  sediment  samples  collected  at  die  three 
locations  mentioned  above,  a  value  of  og  -  1.2  is  used  in  this  study.  In  field  studies,  the 
value  of  og  will  be  estimated  from  a  lognormal  distribution  analysis  of  a  sediment  sample 
from  each  location  of  interest 

To  estimate  the  mean  geometric  particle  radius  of  the  sediment,  die  measured  ratios 
calculated  from  die  right-hand  side  of  Equation  (2.10)  are  compared  with  the  theoretical 
mdosof  F\xi,  og)IF(Xx,og).  In  practice,  the  estimate  of  ag  is  obtained  from  a  match  to 
the  value  of  the  F{Xy  ag)/F(Xl ,  ag)  function  which  is  stored  in  a  lookup  table  in  the 
analysis  software. 

The  mass  concentration  in  the  same  range  bin  can  be  estimated  by  explicitly  solving 
Equation  (2.9)  for  M, 


M 


Ps 


K2 


S?F2(*,,Oj)  SlF2iXy  Og) 


(2.11) 
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Figure  2.6  Theoretical  values  of  F\Xy  o  )lF(Xt ,  a  )as  a  function  of  a,  for  og  =  1.1, 1.2, 
1.3, 1.4, 1.5. 
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where  p,  is  the  density  of  the  scatterers,  K^.  is  the  digital  backscatter  counts  for  each 
frequency,  S  is  die  sensitivity  for  each  frequency,  and  F(X,  is  calculated  from  Equation 
(2.10). 

This  method  therefore  provides  a  technique  for  estimating  mass  concentration  and  the 
mean  geometric  particle  radius  from  the  dual  frequency  backscatter  levels  and  the  assumed 
lognormal  sediment  distribution  standard  deviation  parameter  ag.  A  method  for  determining 
the  system  sensitivity,  S,  for  each  frequency  of  the  CASP  is  discussed  next 

D.  SYSTEM  SENSITIVITY 

Following  the  theoretical  development  of  Sheng  and  Hay  (1993),  the  estimate  of  the 
sensitivity  for  a  frequency  using  wires,  Sw  can  be  expressed  as 


(2.12) 


where  r  is  die  perpendicular  distance  from  the  transducer  to  die  wire,  1/  r  is  die  radial 
spreading  term,  are  the  digital  counts  of  the  backscatter  amplitudes,  and  k  is  die 
compression  wave  number.  The  index  of  summation,  i,  is  determined  by  the  number  of 
different  sized  standard  targets  used.  Four  stainless  steel  wires  of  different  diameters  were 
used  in  this  study. 
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For  a  pulsed  wave  acoustic  source,  the  function  r,  describes  die  dependence  of  the 
backscatter  amplitude  on  the  effective  length  of  the  wire  and  is  expressed  by  (Sheng  and 
Hay,  1993) 


(2.13) 


where  (3,  and  0,  are  functions  of  rt ,  the  distance  to  the  leading  edge  of  the  acoustic  wave 
front,  and  r0  is  the  perpendicular  distance  from  the  transducer  to  the  fixed  wire. 


p,-0f  P2=arocoa(r0/r1),  r#*r,sr0«e cp., 

P,  =  0.  P2  =  P,..  ^secp.sr^ro  +  Wct 


A  schematic  of  the  coordinate  system  for  the  pulsed  system  is  shown  in  I-igure  2.7.  Figure 
2.8  shows  r  as  a  function  of  the  parameter  the  ratio  of  the  radius  of  the  main  lobe  of 

die  transducer  directivity  pattern,  r0tan  ,  to  the  radius  of  die  first  Fresnel  zone 


♦  (r#)  = 


r0t»P- 


t«p. 


(2.14) 


N  2r°x 


where  t  depends  on  the  acoustic  frequency  and  beam  width.  The  values  of  *(/•„)  and  r, 
for  the  1.3  MHz  and  S.2  MHz  transducers  in  this  study  are  0.98,  0.S;  and  2.2,  0.7 
respectively. 


I 

; 
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Figure  2.7  Schematic  of  the  coordinate  system  used  for  an  acoustic  pulse  incident  on  a  fixed 
wire  of  infinite  length.  T  represents  the  transducer.  The  other  symbols  are  explained  in  the 
text 


22 


I  ggjnma  (3 


08 
0.6 
4 
2 
0 

0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5  5.5  6 

psi 


Figure  2.8  IJ,  as  a  function  of  W.rJ  for  t  monostatic  pulsed  acoustic  system. 
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The  form  factor  for  an  infinite  cylinder,  /„  (x),  can  be  expressed  by  (Sheng  and  Hay, 


1993) 

00 

f"  £  (-rf€B«n%«q>(  -^ll)  (2.15) 

ifix 

where  the  Neumann  factor  ^  =  1  forn~Oand  €n  =  2  for  »  *  1  ,  and  is  the  phase  shift 
of  the  nth  partial  wave  (Faran,  1951).  Figure  2.9  shows  the  dependence  of  the  form  factor 
on  the  acoustic  frequency  and  the  radius  of  the  wire.  Table  2  lists  the  values  of  kug  and  /-<*) 
ftn- this  study. 


TABLE  2  VALUES  OF  fa,  AND  /Jr)  FOR  EACH  FREQUENCY  AND  WIRE  RADIUS 


Wire 

1.323 

MHz 

5.292 

MHz 

Radius 

(microns) 

^  On'') 

/JM 

(in') 

/ooW 

25 

■391 

0.12 

0.56 

■SB 

40 

0.22 

0.22 

0.90 

0.92 

50 

0.28 

0.32 

0.12 

0.81 

64 

0.36 

0.44 

1.44 

0.72 

The  relationship  of  /„(x)  to  ka  for  1.3  MHz  is  generally  linear  while  die  relationship  for 
5.2  MHz  shows  the  values  of  /„(x)  increase  for  the  first  two  values  of  ka  then  decrease  for 
die  remaining  two  value.  This  finequency/target  size  dependency  of  /.(x)  allows  the  use  of 
a  multi-frequency  approach  to  discriminate  between  different  sized  targets. 
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IIX  LABORATORY  ACOUSTIC  BACKSCATTER  MEASUREMENTS 


The  absolute  backscatter  sensitivity  for  each  frequency  of  the  CASP  was  calculated  uang 
die  procedure  discussed  in  Section  II D.  The  standard  targets,  four  stainless  steel  wires  of 
different  diameters,  wore  suspended  in  a  600  gallon  tank  while  the  primary  sample  volume 
of  the  CASP  was  moved  through  a  grid  of  points  centered  on  the  wire.  The  CASP  was 
paused  at  each  ghd  point  and  backscatter  levels  wot  measured  to  calculate  the  sensitivity 
for  each  frequency.  The  backscatter  levels  were  also  used  to  compare  the  monostatic  beam 
characteristics  of  the  width  of  die  find  null  point,  half  power  bandwidth,  and  pulse 

length  to  theoretical  values. 

The  inversion  technique  discussed  in  die  previous  chapter  was  verified  by  developing  a 
method  to  measure  die  backscatter  levels  of  known  sediment  concentrations  in  a  unique  test 
voaseL  Clean  and  sorted  sediment  samples  of  increasing  cumulative  weight  were  placed  into 
a  known  volume  of  water  in  a  test  vesseL  The  sediment  was  suspended  by  continuous 
stirring,  then  ensonified  by  die  CASP  to  measure  the  backscatter  levels  from  die  increasing 
concentrations. 

A.  COHERENT  ACOUSTIC  SEDIMENT  PROBE  (CASP) 

The  CASP  is  a  robustly  constructed,  multi-frequency,  underwater  instrument  package 
originally  designed  to  conduct  studies  of  turbulence  in  die  surface  boundary  borer.  It  has 
three  5.2  MHz  acoustic  transducers  placed  radially  around  a  single  1.3  MHz  transducer,  a 
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two  axis  tilt  sensor,  and  a  three  axis  accelerometer  to  detect  the  instrument's  orientation  and 
motion.  In  the  biatatic  mode,  one  of  the  three  5.2  MHz  transducers  emits  a  series  of  32  short 
(O.S  cm)  pulses  into  the  water  column,  while  the  other  two  transducers  receive  the 
backacattered  acoustic  energy  from  die  sample  volume  defined  by  the  intersection  of  the 
transmitter  and  receiver  beams.  After  32  backacattered  pulses  have  been  sampled  by  both 
receiving  transducers,  the  transducers  are  electronically  switched  so  another  S.2  MHz 
transducer  emits  a  pulse  and  the  other  two  receive.  The  monostatic  sampling  mode  is 
interleaved  into  the  bistatic  sampling  sequence  to  estimate  the  along-beam  backscatter 
intensity  for  die  three  S.2  MHz  transducers  and  die  single  1.3  MHz  transducer.  In  die 
monostatic  mode,  a  triplet  of  longer  (1.68  on)  pulses  are  transmitted  from  each  transducer 
in  sequence.  As  each  pulse  is  transmitted,  die  transducer  is  switched  to  the  receive  mode  and 
die  backscatter  amplitudes  are  sampled  in  each  1.68  cm  range  bin  out  to  1.2  meters.  The 
complete  sequence  of  three  bistatic  and  three  monostatic  transmit  modes  is  continuously 
repeated  at  a  rate  of  35  Hz  (Stanton,  1993). 

The  1.3  MHz  transducer  is  mounted  flush  in  die  center  of  die  CASP  instrument  head,  and 
the  5.2  MHz  transducers  are  mounted  cm  10  cm  arms  arranged  radially  outward  from  die 
instrument  head  at  0°,  120°,  and  240°  (where  0°  denotes  die  top  of  die  instrument  head). 
Figprc  3.1  illustrates  the  acoustic  beam  geometry  of  both  die  bistatic  and  monostatic  inodes. 
The  four  acoustic  beams  intersect  at  a  primary  sample  volume  25.4  cm  downbeam  from  die 
1.3  MHz  transducer. 
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Figure  3. 1  Schematic  of  the  CASP  acoustic  beam  geometiy  for  the  bistatic  and  monos talic 
modes.  (Stanton  1993) 
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B.  ABSOLUTE  HACKSCATTER  SENSITIVITY 

The  absolute  system  reaponac  of  the  CASP  to  backscatlcr  from  cylindrical  targets  of  known 
dwmeter  waa  calculated  following  the  procedure  diacuwod  in  Section  n  D.  Previous  studies 
have  been  based  on  spherical  targets,  particularly  lead-glass  beads.  Spheres  have  the 
advantage  that  their  target  strength  is  independent  of  orientation,  but  it  can  be  difficult  to 
suspend  a  spherical  target  precisely  in  the  center  of  die  acoustic  beam,  particularly  in  die 
very  narrow  acoustic  beams  generated  by  transducers  operating  at  MHz  frequencies.  Woe 
targets  offer  die  advantage  that  the  mounting  structure  is  located  completely  outside  the 
acoustic  beam  path  so  backscatter  from  the  supporting  structure  is  not  present 

Sheng  and  Hay  (1993)  conducted  a  study  comparing  system  sensitivity  calculations  usn% 
lead-glass  beads  with  those  calculated  using  stainless  steel  wires.  There  was  good  agreement 
between  die  two  approaches  when  an  appropriate  correction  factor  was  made,  indicating  that 
stainless  stcd  wires  can  be  used  as  standard  targets  in  system  sensitivity  calculations. 

In  this  experiment,  system  sensitivity  was  determined  using  stainless  steel  wires  of  four 
different  retfii:  25, 40,  SO,  and  64  microns.  Table  3  fists  die  values  of  x  (=  kpS )  for  c  =  1480 
m/s  for  each  wire  and  frequency  combination.  The  values  of  k,  for  die  1.3  MHz  and  S.2 
MHz  transducers  were  S617  m1  and  22467  m1  respectively. 
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TABLE  3  RAM  AND  VALUES  OF  x  =  k*(mJ)  FOR  EACH  FREQUENCY  AND  WIRE 
COMBINATION 


Radius  (pM) 

1.323  MHz 

5.292  MHz 

25 

0.14 

0.36 

40 

0.22 

30 

0.28 

1.12 

64 

0.36 

1.44 

Each  wire  was  suspended  under  tension  in  the  test  tank  The  wires  were  oriented  parallel 
to  the  transducer  face  which  was  being  investigated  (vertical  for  1.3  MHz,  23.$°  from  die 
vertical  for  3.2  MHz).  The  CASP  was  moved  under  computer  control  in  small  increments 
in  X/Y  space  by  a  digital  servo  controller  which  positioned  a  precision  stage  supporting  die 
CASP.  The  servo  controller  is  capable  of  incremental  movements  of  the  CASP  staging  as 
small  as  30  microns.  The  focal  point  of  die  acoustic  beans  traced  a  grid  of  points  with  the 
wire  located  approximately  in  die  grid  center.  Backscatter  levels  woe  measured  by  die 
CASP  at  each  grid  point  for  a  specified  period  of  time.  Figure  3.2  shows  a  plan  view  of  die 
gird  of  points  used  in  die  absolute  backscatter  runs.  The  grid  dimensions  in  these 
calibrations  were  3  cm  by  10  cm  with  an  increment  of  0.2  on.  The  dwell  time  at  each  grid 
point  was  10  seconds. 
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Figure  3.2  Plan  view  of  the  5  cm  by  10  cm  grid  traced  out  by  the  primary  sampling  volume 
of  the  CASP  during  the  absolute  backscatter  runs.  The  data  collection  nm  began  with  the 
primary  sampling  volume  in  the  (0,  0)  position  with  the  arrows  indicating  the  motion  of  the 
CASP.  The  solid  arrows  indicate  data  collection,  while  the  dashed  arrows  indicate 
repositioning  of  the  CASP  in  increments  of  0.2  cm. 
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C.  LABORATORY  SEDIMENT  BACKSCATTER  TESTS 
1.  Experimental  Tank  and  Mixing  Vessel 
A  600  gallon  acrylic  tank  was  used  to  contain  the  CASP  in  fresh,  filtered  water  during 
the  calibration  procedures.  The  tank  was  equipped  with  a  computer-controlled  X/Y 
positioning  stage  under  which  the  CASP  was  suspended.  The  relative  size  of  die  tank  and 
die  CASP  (Figure  3.3)  allowed  target*  to  be  positioned  in  the  path  of  die  acoustic  beams  as 
the  CASP  stage  was  repositioned.  An  acrylic  test  vessel  was  constructed  with  a  plastic  film 
window  to  minimize  attenuation  across  the  vessel-water  boundary.  A  controllable  speed, 
vertically  oriented  propeller  stirrer  was  positioned  in  a  vertical  acrylic  cylinder  serving  as 
a  directional  nozzle  which  was  secured  to  die  tank  frame  to  center  it  in  the  mixing  vessel. 
The  mixing  vessel  and  propeller  stirrer  nozzle  were  designed  so  that  no  level  surfaces  or 
flow  stagnation  points  were  present  which  could  allow  sediment  to  settle  out  of  die  mixing 
volume  and  reduce  die  test  sample  concentration.  The  mixing  vessel,  containing  15.4  liters 
of  dean,  deaerated  water,  was  placed  in  die  test  tank  with  an  OBS  inserted  in  an  O  ring 
sealed  sidewall  opening  with  die  transducer  face  positioned  downward.  The  temperature  of 
the  tank  and  mixing  vessel  water  were  measured  to  0.1°  C  accuracy.  The  tank  and  vessel 
temperatures  dfflfered  by  at  most  1°  C  throughout  the  series  of  test  runs.  Any  bubbles  on  the 
CASP  or  mixing  vessel  surfaces  were  brushed  off,  with  particular  care  being  taken  around 
the  transducer  faces  and  mixing  vessel  window.  The  CASP  was  positioned  such  that  the 
acoustic  beam  intersection  point  was  approximately  1.0  cm  inside  the  mixer  vessel  window. 

The  propeller  stirrer  nozzle  was  placed  downward  within  one  degree  of  vertical  in  die  test 
vessel  and  secured  to  the  tank  frame.  The  calibration  run  for  each  grain  size  range  was 
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1.  COHERENT  ACOUSTIC  SEDIMENT  PROBE  (CASP) 

2.  ACOUSTICALLY  TRANSPARENT  MEMBRANE 

3.  ACRYLIC  NOZZLE  AM)  DUCTED  PROPELLER 

4.  ACRYLIC  TEST  VESSEL 

5.  VARIABLE  SPEED  MOTOR 

«.  OPTICAL  BACKSCATTER  SENSOR 

7.  ACRYLIC  LABORATORY  ACOUSTIC  TEST  TANK 

8.  SUSPENDED  SEDIMENT 


Figure  3.3  Schematic  of  the  test  tank,  CASP,  and  the  mixing  vessel  with  propeller  stirTer. 
(Stanton,  1993) 


performed  by  adding  sediment  samples  to  the  test  vessel  with  the  propeller  stirrer  running, 
resulting  in  a  homogeneous  mass  concentration.  An  equal  amount  of  clear  water  from  the 
vessel  was  removed  in  the  sample  vial  to  maintain  the  mass  concentration  for  that  particular 
rua  Backscatter  data  were  collected  for  120  seconds.  The  sorted  samples  woe  added  to  the 
vessel  in  sequence,  with  the  final  weight  of  sediment  in  die  test  vessel  totaling  ?2  grams. 
This  gave  a  maximum  mass  concentration  of  2.08  g/1. 


2.  Sediment  Collection  and  Preparation 

Sediment  samples  woe  obtained  from  three  locations  along  the  Monterey  Bay  coast: 
Arikmar  State  Beach,  Naval  Postgraduate  School  Beach,  and  die  Monterey  Wharf.  Samples 
from  each  location  were  cleaned,  sorted,  and  weighed  to  0.01  gram  accuracy.  Sieves  with 
mesh  sizes  of 0.090  mm,  0.125  mm,  0.180  ram,  0.250  mm,  0.355  mm,  1.00  mm,  and  2.00 
mm  woe  used  in  die  sorting  of  the  samples.  Discrete  samples  of  clean  and  sorted  sediment 
wwe  weighed  into  vials  containing  a  weak  solution  of  water  and  Triton  100.  Sorted  samples 
woe  weighed  out  such  that  sequential  additions  of  the  sediment  samples  into  the  test  vessel 
readied in cumdative suboctave  weights  of  0.5, 1.0, 1.5,  2.0, 4.0,  6.0,  8.0, 12.0, 16.0, 24.0, 
and  3Z0  grams.  Triton  100  is  a  wetting  agent  that  prevented  the  sediment  grains  from  being 
encased  in  air  bubbles  wh»oh  could  significantly  increase  die  backscatter  cross  section 
resulting  in  erroneous  backscatter  levels. 

Backscatter  data  wen  analyzed  from  the  first  range  bin  inside  die  mixing  vessel 
window.  Figures  3.4,  3.5,  and  3.6  show  die  mean  monostatic  backscatter  levels  for  each 
range  bin  over  die  duration  of  die  sample  run  (120  sec).  The  left  and  right  panels  show  the 
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<Log(Variance  of  Counts)>  <Log(Vinance  of  Counts)> 


Range  (cm)  Range  (cm) 


Figure  3.4  Mean  monostotic  log  backscatter  power  level  profile  Tor  sediment  of  size  1 25  to 
180  microns.  The  top  panel  shows  mean  backscatter  levels  for  no  sediment  and  no  stirring 
of  the  mixer  vessel  water  volume.  The  bottom  panel  shows  mean  backscatter  levels  for  no 
sediment  with  stirring  of  the  mixer  vessel  water  volume.  The  left  panel  shows  the 
backscatter  level  profile  for  1.3  MHz.  The  right  panel  shows  the  backscatter  level  profile 
for  5.2  MHz. 
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<log(Variance  of  Counts)>  <Log(Variance  of  Counts)> 
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Grain  S\2b:  .0,.  1 25  min. 
Sample  Weigtil:  0.5  g. 
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Figure  3.5  Same  as  Figure  3.4  but  for  0.5  g  and  8.0  g  of  sediment  in  suspension. 
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<Log(Variance  of  Counts)>  <Log(Vanance  of  Counts)> 


Figure  3.6  Same  as  Figure  3.4  but  for  16.0  g  and  32.0  g  of  sediment  in  suspension. 


37 


backscatter  levels  for  the  1.3  MHz  and  5.2  MHz  acoustic  beams  respectively  for  the  grain 
size  diameter  range  of  125  to  180  microns.  The  cumulative  sample  weight  of  sediment 
shown  am  0  grams  with  no  mixing  of  the  water  volume,  and  0,  0.5,  8.0,  16.0  and  32.0  grams 
widt  mixing  of  die  water  column.  For  both  the  1.3  and  5.2  MHz  frequencies,  the  0  gram,  no 
mixing  case  exhibits  backscatter  levels  which  can  be  attributed  to  micro-turbulent  vortices 
and  small  contaminants  present  in  the  test  tank  and  mixer  water  volumes  which  represent  a 
noise  floor  for  this  calibration  method.  For  the  mixing  cases  the  5.2  MHz  profiles  show  a 
kage  increase  in  backscatter  power  (by  a  factor  of  10s  for  32  grams)  in  the  mean  backscatter 
power  returns  inside  die  mixer  window  (starting  at  26  cm  range)  at  each  concentration.  The 
mean  backscatter  returns  for  the  1.3  MHz  beam  also  increase  inside  the  mixer  window  (by 
a  factor  of  neariy  10*  for  32  grams)  for  die  higher  concentrations  but  continue  to  increase 
for  the  lower  concentrations,  to  reach  a  maximum  at  a  range  of  about  32  cm.  The 
backscatter  drops  off  sharply  at  about  34  cm  where  the  backscatter  profile  resembles  more 
closely  that  of  die  5.2  MHz  acoustic  beam.  The  higher  values  of  backscatter  levels  for  die 
low  concentration  profiles  for  1.3  MHz  near  die  20-25  cm  range  are  believed  to  be  due  to 
the  plastic  film  window.  Backscatter  values  are  sampled  in  the  range  bins  through  die 
suspended  sediment  volume. 

Appendix  A  contains  the  sieve  analysts  for  each  of  die  nine  samples  using  die  measured 
weights  for  each  sieve  size  range.  The  analyses  show  die  size  distributions  were 
significantly  different  for  each  of  the  three  locations.  The  mean  percentages  of  sediment  size 
ranges  for  die  three  collection  sites  shown  in  Table  4  indicate  NPS  Beach  Lab  sediment 
consists  mainly  of  medium  and  fine  sized  silts  with  most  sediment  between  0.063  and 
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0.5  nun.  Monterey  Wharf  sediment  is  similar  to  die  sediment  collected  at  Asilomar  but  has 
a  more  even  proportion  of  larger  sued  grains  in  the  0.5  to  1.0  mm  range. 


TABLF  J  PERCENTAGE  WEIGHT  BY  SEDIMENT  CLASSIFICATION 


Classification 

Size  Range  (mm) 

Asilomar 

Beach  Lab 

M.  Wharf 

Gravel 

>1.0 

17.02 

19.35 

Coarse  Sand 

25.81 

1.26 

15.04 

Med./Fme  Silt 

0.063  -  0.500 

56.15 

97.16 

64.09 

Coarse-Fine  Silt 

<0.063 

1.02 

1.58 

1.52 
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IV.  RESULTS  AND  DISCUSSION 


The  objectives  pursued  in  this  study  were  the  determination  of  the  sensitivity  of  each 
frequency  of  the  CASP,  and  development  and  verification  of  a  two-frequency  algorithm 
inverting  backscattered  acoustic  energy  from  a  volume  of  suspended  sediment  to  estimate 
mass  concentration  and  a  geometric  mean  radius  of  the  suspended  sediment.  The  inversion 
assumes  the  sediment  is  described  by  a  lognormal  size  distribution  with  a  prescribed 
variance  aftiie  distribution  about  the  mean  radius.  The  system  sensitivity  for  each  frequency 
was  determined  from  tire  backscatter  levels  of  four  stainless  steel  wires  of  known  diameters 
using  the  procedure  outlined  in  Section  III  B.  Directivity  pattern  characteristics  obtained 
from  graphical  analysis  of  the  backscatter  levels  from  the  four  wires  were  compared  to  the 
theoretical  values  of  rite  monostatic  mode  for  each  CASP  frequency. 

A.  WIRE  TARGET  SYSTEM  SENSITIVITY  ESTIMATION 

The  system  sensitivity  for  the  CASP  was  determined  by  ensorafying  each  wire  separately, 
following  the  procedure  discussed  in  Section  HI  B.  Table  S  lists  die  rms  counts,  tire  values 
of  ka,  Ae  form  factor,  the  sensitivity  for  each  wire;  and  the  overall  sensitivity  calculated  for 
the  1.3  MHz  transducer.  Table  6  hats  the  same  information  for  the  S.2  MHz  transducer. 
Pram  Table  S  and  Table  6,  the  values  of  ka  ranges  from  0.14  to  0.36  for  1.3  MHz  and  0.56 
to  1.44  for  5.2  MHz.  The  form  factor  has  the  range  of  values  of  x  from  0.14  to  0.36  and  is 
generally  linear,  as  is  the  plot  of  peak  count  values  for  1.3  MHz.  The  form  factor  for  the 
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TABLE  5  PARAMETERS  FOR  1.323  MHZ  SENSITIVITY  CALCULATIONS 


1.323MHz 

Wire  Radius 
(microns) 

25 

40 

50 

64 

RMS  Counts 

1694*92 

6194*174 

19929*1066 

30500*2770 

l 

I! 

w 

X 

0.14 

0.22 

0.28 

0.36 

LCttl 

0.12 

0.22 

0.32 

0.44 

Numerator  (counts) 

mm 

5.75 

30.1 

71.7 

Denominator  (10"6) 

0.16 

.86 

2.28 

5.52 

S  (10*  counts) 

4.23*0.30 

6.67*0.33 

13.2*0.92 

13.0*1.43 

<N> 

108.2 

<D>(10*) 

8.8 

Sw  (10^  counts) 

9.28*5.28 

TABLE  6  PARAMETERS  FOR  5.292  MHZ  SENSITIVITY  CALCULATIONS 


5.292MHz 

Wire  Radius 
(microns) 

25 

40 

50 

64 

RMS  Counts 

10922*884 

20557*3575 

17108*1444 

16395*1186 

Xj(  -ka) 

0.56 

0.90 

1.12 

1.44 

IX  60  f 

0.74 

0.92 

0.81 

0.72 

Numerator 

37.77 

111.78 

91.58 

88.26 

Denommator  (10*) 

11.96 

29.57 

28.65 

28.65 

S  (10*  counts) 

3.16*0.28 

3.78*0.68 

3.20*0.29 

3.05*0.24 

<N> 

329.4 

0(10*) 

99.2 

Sw  (10*  counts) 

3.30*0.70 
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range  of  vahie  for  5.2  MHz  increases  to  a  relative  maximum  at  a  value  of  x  at  about  0.85, 
then  decreases  to  a  relative  minimum  at  a  vahie  of  x  of  about  1.4. 

The  rats  count  values  for  each  frequency  and  wire  combination  are  plotted  in  Figure  4.1. 
The  sensitivity  values  of  each  wire  and  frequency  combination  are  plotted  in  Figure  4.2. 
Theoretically,  die  sensitivity  for  each  wire  and  a  single  frequency  should  be  the  same  since 
the  form  factor  would  normalize  the  different  response  from  wires  of  different  radii.  The 
plot  of  the  5.2  MHz  sensitivities  for  die  four  wires  shows  this  is  the  case,  with  die 
sensitivities  consistent  to  21%.  The  sensitivities  for  the  1.3  MHz  transducer  however  are 
widely  spread,  suggesting  s  problem  with  die  measurement  procedure.  It  is  suspected  that 
there  is  a  strong  dependence  between  the  angle  made  between  the  wire  and  the  transducer 
face  and  die  backscatter  response,  depute  attempts  to  keep  die  wires  orthogonal  to  die 
transducers.  The  wide  range  of  values  for  die  1.3  MHz  transducer  prevented  further  use  of 
die  wire-derived  sensitivity  values.  In  future  work,  a  spherical  target  calibration  wiO  be 
attempted,  although  it  introduces  a  vertical  dependence  to  the  backscatter  values. 

a  MULTI-FREQUENCY  INVERSION  ALGORITHM 

Thirteen  mass  concentrations  for  aedanent  collected  from  Monterey  Wharf  were  used  to 
teat  the  CASP  response  to  a  known  sediment  size  distribution  and  controlled  concentrations. 
Figures  4.3  shows  the  mean  monostadc  backscatter  level  profile  for  each  range  bin  averaged 
over  die  duration  of  die  run  for  a  mass  concentration  of  0.26  g/L  The  top  panel  shows  die 
profife  for  the  1.3  MHz  transducer  while  the  bottom  panel  shows  the  profile  for  the  5.2  MHz 
transducer.  Figure  4.4  shows  die  profile  fra*  a  mass  concentration  of  2. 1  g/1.  Regression 
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Sensitivity  (counts) 


x  106 


Figure  4.2  Sensitivity  values  for  each  wire  and  frequency  combination 
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Figure  4.3  Monostotic  backscatter  level  profile  for  each  range  bin  for  a  mass  concentration 
of  0.26  g/1.  The  top  and  bottom  panels  show  the  profiles  for  the  1.3  MHz  and  5.2  MHz 
transducers  respectively. 


5 


Counts 


Figure  4.4  Same  as  Figure  4.3  but  for  a  mass  concentration  of  2. 1  g/1. 
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knot  of  the  backscatter  levds  for  bin  13, 14,  and  IS  for  die  1.3  MHz  and  5.2  MHz  transducer 
have  been  drawn  in.  The  backacatter  level  of  range  bin  12,  the  first  range  bin  inside  the 
mixing  vessel,  fits  the  regression  line  generally  well  for  the  S.2  MHz  transducer  while  the 
ban  12  backscatter  level  for  the  1.3  MHz  transducer  is  seen  to  be  well  above  the  extrapolated 
value  from  the  regression  line. 

As  previously  discussed  in  Section  1C2,  these  largo-  then  expected  backacatter  levels 
may  be  attributed  to  scatter  by  the  plastic  film  window  or  multiple  backacatter  from  the 
transducer  beam  side  lobes  coming  into  contact  wife  the  mixing  vessel  window  supporting 
structure.  To  minimize  the  influence  of  this  multiple  backscatter  from  the  mixing  vessel  on 
the  analysis  of  die  data,  the  bin  12  backacatter  level  for  the  1.3  MHz  transducer  bin  was 
extrapolated  from  the  regression  for  each  run.  The  bin  12  backacatter  level  for  the  0  gram, 
mixing  run  was  considered  the  baseline  backscatter  value  and  was  subtracted  from  the 
backacatter  level  for  each  subsequent  run .  The  bin  12  backscatter  levels  for  the  5.2  MHz 
transducer  did  not  require  correction  and  were  analyzed  as  measured. 

Table  7  lists  die  cumulative  mass,  cumulative  mass  concentration,  die  data  archive  file, 
and  the  mean  digital  counts  of  die  backacatter  amplitudes  for  both  frequencies.  The 
backacatter  levels  for  the  1.3  MHz  transducer  have  been  corrected  as  described  above. 


\ 

i 
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TABLE  7  BROAD  SPECTRUM  SUSPENDED  SEDIMENT  DATA 


Cum.  Mass 
(2) 

Cum.  Cone. 

m 

CASPTS 

1.3  MHz 
(counts) 

5.2  MHz 
(counts) 

0 

0 

55 

-29*0 

-151*26 

0 

0 

56 

0 

0 

0.5 

0.03 

57 

59*18 

1704*204 

1.0 

0.65 

62 

979*39 

2409*265 

1.5 

0.10 

66 

1307*54 

2991*239 

2.0 

0.13 

67 

1613*58 

3534*318 

3.0 

■331 

68 

2066*99 

4293*343 

4.0 

0.26 

69 

2604*120 

5034*403 

6.0 

0.39 

58 

3270*177 

6279*565 

8.0 

0.52 

59 

3968*238 

7329*660 

12.0 

0.78 

60 

5109*511 

9139*731 

16.0 

1.0 

61 

9753*858 

10294*926 

24.0 

1.6 

63 

7522*617 

13000*1170 

32.0 

2.1 

64 

8737*664 

14712*1030 

Figure  4.5  shows  the  {dot  of  the  dependence  of  the  backscatter  levels  of  die  suspended 
sedtuent  from  Monterey  Wharf  to  the  cumulative  concentration.  From  Equation  (2.9),  the 
backacatter  levels  are  proportional  to  the  square  root  of  the  mass  concentration.  One  would 
expect  a  linear  relationship  to  arise  from  die  square  of  the  backscatter  levels  as  a  function 
of  die  mass  concentration.  Figure  4.6  shows  the  plot  of  die  square  of  the  backscatter  levels 
for  each  cumulative  concentration.  Ignoring  die  obvious  outfcer  of  the  1.3  MHz  backscatter 
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Figure  4.5  Plot  of  the  backscatter  levels  as  a  function  of  mass  concentration. 
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Figure  4.6  Plot  of  the  square  of  the  backscatter  levels  as  a  function  of  mass  concentration. 
The  linear  nature  of  the  form  verifies  the  square  root  proportionality  between  backscatter  and 
mass  concentration. 
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forte  1.0  gl  mass  concentration,  believed  to  be  caused  by  a  system  malfunction,  Figure  4.6 
does  indeed  verify  this  linear  relationship  within  the  standard  deviations  of  the  backscatter 
counts  for  each  mass  concentration.  This  satisfying  result  emphasizes  foe  reliability  of  foe 
equipment  and  procedures  toed  in  this  part  of  foe  study. 

The  slope  of  foe  regression  lines  discussed  above  represents  foe  decrease  in  backscatter 
level  with  distance,  or  attenuation,  for  each  sample  concentration,  hi  this  case,  it  specifically 
represents  foe  attenuation  ova*  the  length  of  one  range  bin  (1.68  cm)  far  foe  CASP.  Table 
8  lists  the  attenuation  for  each  cumulative  concentration  and  frequency.  Figure  4.7  plots 
•ttmutinn  versus  concentration  values  which  can  be  used  in  field  measurements  to  estimate 
the  attenuation  along  each  beam  in  foe  presence  of  sediment  Starting  at  range  bin  1,  foe 
mass  concentration  can  be  estimated  using  Equation  (2.1  IX  foot  the  local  attenuation  can 
be  estimated  from  a  lookup  table  of  foe  form  of  Table  8  (again,  foe  attenuation  for  the  1.3 
MHz,  1  gdmass  concentration  should  be  ignored).  The  backscatter  level  for  foe  subsequent 
bins  can  therefore  have  foe  effects  of  attenuation  removed. 
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Figure  4.7  Attenuation  as  a  function  of  mass  concentration. 


TABLE  8  ATTENUATION  FOR  EACH  CUMULATIVE  CONCENTRATION 


Mass  Concentration 

Cg/I) 

1.3  MHz 
(103  counts/m) 

5.2  MHz 
(lO’counts/m) 

0 

-4.38 

-.00324 

0.03 

-2.97 

-4.76 

0.65 

-3.41 

-6.71 

0.10 

•4.06 

-8.24 

0.13 

-4.62 

-9.97 

0.19 

-4.71 

-12.1 

0.26 

-6.91 

-15.3 

0.39 

-7.91 

-19.7 

0.52 

-10.1 

-25.8 

0.78 

-13.8 

-33.5 

1.0 

-40.4 

-40.41 

1.6 

-20.5 

-53.7 

2.1 

-24.6 

-67.5 

As  the  wire  sensitivity  calibrations  hive  unresolved  problems  for  the  1.3  MHz  transducer, 
the  direct  measurements  of  backscatter  versus  a  known  sediment  concentration  shown  in 
Table  7  and  Figure  4.S  can  be  used  to  estimate  die  sensitivities  S,  and  S}.  Following 
Equation  (2.11),  the  sensitivity  can  be  expressed  as 

s  =  F  *•'  <41> 
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for  each  concentration,  j,  as  aO  these  variables  are  known.  For  concentrations  well  above 
the  cattxabon  tank  noise,  and  by  rejecting  the  clear  outlier  for  the  1.3  MHz  transducer  at  die 
mass  concentration  of  1.0  g/1,  the  sensitivities  in  Table  9  are  obtained. 


TABLE  9  SENSITIVITIES  FOR  EACH  C ASP  TRANSDUCER  CALCULATED  FROM 
MONTEREY  WHARF  SEDIMENT  PARAMETERS 


Cumulative 

Concentration 

Sensitivity 

(counts) 

(g/1) 

1.3  MHz 

5.2  MHz  A 

5.2  MHz  B 

5.2  MHz  C 

0.19 

15208 

16112 

14720 

16397 

0.26 

15726 

16268 

14538 

16318 

0.39 

15423 

16457 

14440 

16475 

0.S2 

15701 

16571 

15021 

16430 

0.78 

15948 

16794 

14809 

16646 

1.0 

24840 

16347 

14558 

16312 

1.6 

15956 

16798 

14940 

16348 

2 

15859 

16438 

14585 

16301 

MeaitfcSTD 

15867±591 

16421±237 

14698±177 

16475*218 

Fjpve  4.8  shows  the  plot  of  the  values  of  die  sensitivities,  calculated  using  Equation  (4.1) 
and  Ihe  sediment  {Nnmeters  for  Monterey  Wharf  for  die  1.3  MHz  and  each  of  the  three  $.2 
MHz  transducers.  The  calibration  tank  noise  is  seen  in  die  apparent  high  sensitivity  values 
for  suspended  sedhnent  concentrations  less  than  0.19  g/1.  The  sensitivity  values  for  the 
concentrations  equal  or  greater  than  0.19  g/1  indicate  a  high  consistency  throughout  die 
range  of  mass  concentrations  for  each  transducer.  The  measured  sensitivity  range  of  the  1.3 
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MHz  transducer  varies  3.7%,  while  the  sensitivities  for  the  S.2  MHz  transducers  vary  only 
1.4%,  1.2%,  and  Irrespectively. 

Lastly,  the  measured  sample  volume  defined  by  the  transducer  length  and  beam  width  is 
compared  with  the  theoretical  values.  Table  10  lists  the  measured  values  of  fc,,  the  angle 
of  die  first  null  point  of  the  acoustic  beam  directivity  pattern.  The  theoretical  values 
calculated  in  Equation  (2.1)  were  3.1°  for  1.323  MHz  and  3.1°  for  3.292  MHz. 


TABLE  10  ANGLE  OF  THE  FIRST  NULL  POINT 


Frequency 

(MHz) 

IP«1  (degrees) 

Wire  Radius  (microns) 

<Pml> 

25 

40 

50 

64 

1.323 

2.04 

2.66 

2.66 

2.66 

2.51 

5.292 

2.86 

2.86 

2.86 

2.65 

2.81 

Figures  4.9  to  4.24  show  3  D  representations  of  the  directivity  pattern  for  the  1.3  MHz 
and  5.2  MHz  monostatic  acoustic  pulse  for  bin  12.  The  X  and  Y  plane  represei  j  the  grid 
the  CASP  was  moved  through  while  badkscatter  data  were  collected.  The  Z  axis  represents 
the  counts  and  die  log  of  die  square  of  die  counts  (as  generated  by  die  CASP  processing 
program)  as  labeled.  Also  shown  are  end  and  side  views  of  the  acoustic  pulse  in  log  space. 

Table  11  lists  the  graphically  measured  values  of  the  half-power  bandwidth.  The 
theoretical  value  for  both  1.3  and  5.2  MHz  was  2.6°. 
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Figure  4.9  3  D  representation  of  the  monostatic  acoustic  pulse  for  1.3  MHz  and  wire  of 
radius  25  microns  in  log  and  linear  space. 
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Figure  4. 10  End  and  side  views  of  the  monostatic  acoustic  pulse  for  1 .3  MI  1/  and  wire  of 
radius  25  microns. 


Figure  4.1 1  3D  representation  of  the  monostatic  acoustic  pulse  for  1.3  MHz  and  wire  of 
radius  40  microns  in  log  and  linear  space. 
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1.3  MHz  Wire  Radius:  40  microns  End  View 
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Figure  4.12  End  and  side  views  of  the  monostatic  acoustic  pulse  for  1.3  MHz  and  wire  of 
radius  40  microns. 
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Figure  4.14.  End  and  side  views  of  the  monostatic  acoustic  pulse  for  1.3  MHz  and  wire  of 
radius  SO  microns. 
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Figure  4.15  3  D  representation  of  the  monostatic  acoustic  pulse  for  1.3  MHz  and  wire  of 
radius  64  microns  in  log  and  linear  space. 


Y  Grid  (cm) 


Figure  4. 16.  End  and  side  views  of  the  monostatic  acoustic  pulse  for  1 .3  Ml  Iz  and  wire  of 
radius  64  microns. 
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Figure  4.17  3D  representation  of  the  monostatic  acoustic  pulse  for  5.2  MHz  and  wire  of 
radius  25  microns  in  log  and  linear  space. 
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Figure  4.18  End  and  side  views  of  the  monostatic  acoustic  pulse  for  S.2  MHz  and  wire  of 
radius  25  microns. 
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Figure  4.19  3  D  representation  of  the  monostatic  acoustic  pulse  for  5.2  MIIz  and  wire  of 
radius  40  microns  in  log  and  linear  space. 
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Figure  4.20  End  and  side  views  of  die  monostatic  acoustic  p»  or  5.2  MTf?  and  wire  of 
radius  40  microns. 
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Figure  4.21  3  D  representation  of  the  monostatic  acoustic  pulse  for  5.2  MHz.  and  wire  of 
radius  50  microns  in  log  and  linear  space. 
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Figure  4.23  3  D  representation  of  the  monostatic  acoustic  pulse  for  5.2  MHz  and  wire  of 
radius  64  microns  in  log  and  linear  space. 
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Figure  4.24  End  and  side  views  of  die  monostatic  acoustic  pulse  for  $.2  MHz  and  wire  of 
radius  64  microns. 
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TABLE  1 1  HALF-POWER  BANDWIDTH 


Frequency 

(MHz) 

|  |  (degrees) 

Wire  Radius  (microns) 

<ip„i> 

25 

40 

50 

64 

1.323 

1.2 

1.2 

1.2 

1.2 

1.2 

5.292 

1.2 

1.2 

1.2 

1.2 

1.2 

Comparison  of  the  calculated  values  of  ,  and  indicate  the  acoustic  beam  for  both 
1.3  MHz  and  3.2  MHz  are  slightly  narrower  than  calculated  by  theory,  but  show  a  close  to 
ideal  beam  response. 

Table  12  lists  L>  ,  die  length  of  the  monostatic  pulse  length  for  both  frequencies.  The 
theoretical  value  is  0.0168  meten.  Again,  there  is  good  agreement  between  die  calculated 
and  measured  values  suggesting  that  die  acoustic  transmitter,  transducer  and  receiver  system 
is  performing  as  designed . 


TABLE  12  MONOSTAHC  PULSE  LENGTH 


Frequency 

(MHz) 

( )  (10**  meter) 

Wire  Radius  (microns) 

(<M> 

25 

40 

50 

64 

1.323 

1.7 

1.7 

1.7 

1.6 

1.7 

5.292 

1.9 

1.7 

1.9 

1.9 

1.9 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Ascfiscuasedin  Section  IV  A  the  sensitivity  for  each  wire  should  be  equal  for  each  of  the 
frequencies  being  evaluated  baaed  on  the  dependence  of  die  form  factor  on  the  frequency  of 
the  incident  energy  and  the  size  of  the  fixed  wire.  The  sensitivities  calculated  for  die  5.2 
MHz  transducer,  3.30±0.70  x  10*  counts,  were  consistent  to  within  21  %.  The  sensitivities 
calculated  for  die  1.3  MHz  transducer  woe  inconsistent  with  a  mean  value  of  9.28±5.28  x 
10*  counts.  This  laige  variation  is  believed  to  be  due  to  a  strong  angular  dependence 
between  the  orientation  of  die  fixed  wire  target  and  die  1.3  MHz  transducer.  While  use  of 
a  fixed  wire  as  a  standard  target  avoids  the  z-axis  dependence  and  die  difficulty  of  locating 
targets  of  minute  size  in  the  narrow  acoustic  beam,  these  advantages  appear  to  be 
outweighed  by  critical  angular  positioning  requirements.  The  sensitivities  of  the  two 
frequencies  are  parameters  required  in  the  estimation  of  both  the  mean  particle  radius  and 
the  mass  concentration  of  die  suspended  sediment  using  the  inversion  method  discussed  in 
this  wok  The  angular  dependency  between  a  fixed  wire  taiget  and  the  1.3  MHz  transducer 
can  be  avoided  by  die  use  of  a  single  small,  spherical  taiget  such  as  a  bead  or  BB,  suspended 
in  die  acoustic  beam.  Although  there  will  be  no  angular  dependence  between  the  spherical 
target  and  the  transducer  due  to  die  symmetry  of  the  taiget,  care  must  be  taken  to  accurately 
locate  die  taiget  in  the  coder  of  the  very  narrow  acoustic  beam. 
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The  values  of  the  aenaiivitiea  obtained  from  Eq.  (4.1)  and  the  known  suspended  sediment 
parameters  were  highly  consistent  for  the  1.3  MHz  and  the  three  3.2  MHz  transducers 
indicating  the  equipment  and  procedure  woe  highly  reliable. 

The  proportionality  between  die  backscatter  levels  and  the  square  root  of  the  mass 
concentration  confirms  the  reliability  of  the  experimental  method  for  the  collection  of 
backscatter  data  from  the  suspended  sediment  The  cause  of  the  single  outlying  point  in  the 
1.3  MHz  data  is  unknown,  but  is  believed  to  be  a  system  malfunction  and  not  indicative  of 
the  technique  or  noise  performance. 

The  form  of  the  curves  for  the  theoretical  ratio  of  F{X9  a }!FQCx,oJ  ,  shown  in  Figure 
2.6,  shows  a  low  sensitivity  to  the  choice  of  op  and  sufficient  slope  to  allow  a  satisfactory 
estimation  of  the  mean  particle  radius  of  die  lognormal  sediment  population  over  a  range 
from  20  to  100  microns  using  the  two-frequency  inversion  method  described  here.  The  slope 
of  die  curves  approaches  zoo  for  radii  greater  than  100  microns  (200  micron  diameter), 
seriously  reducing  the  ability  of  the  nxs ,  og  )IF(Xl ,  og )  ratio  to  discriminate  mean  radius 
sizes  outside  this  range. 

The  reflective  material  of  the  test  tank  results  in  die  constructive  reflection  of  a  portion 
of  the  acoustic  energy  off  the  sides  and  bottom  of  the  tank.  This  creates  regions  of 
ambjgaouB  backscatter  amplitudes  which  can  be  misinterpreted  as  the  backscatter  from  a  test 
target  or  suspended  sediment  under  investigation.  Additional  analysis  of  these  “range” 
ambiguities  is  required  to  identify  their  location  to  avoid  positioning  of  test  targets  or  die 
mixing  vessel  near  these  positions.  This  becomes  particularly  critical  for  low  concentrations 
and  other  weak  targets. 
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Calibration  rum  using  sediment  samples  with  different  lognormal  distribution  parameters 
will  add  to  the  verification  of  the  inversion  method  by  use  of  die  different 


KXitog  )IF(X , , og )  curves  determined  by  the  value  of  o(.  Furthermore,  the  attenuation 

versus  concentration  characteristics  for  different  sediment  types  win  be  applied  in  along- 
hewn  mass  estimates  to  compensate  for  along-beam  attenuation  in  the  backscatter  levels,  and 
could  possibly  provide  additional  discrimination  of  sediment  types. 

The  two  minute  duration  of  the  suspended  sediment  rum  was  chosen  to  allow  die 
collection  of  s  sufficient  amount  of  data  to  ensure  the  effect  of  statistical  noise  was 
nmmized.  Further  investigation  of  the  temporal  resolution  of  this  inversion  technique  will 
reveal  the  minimum  sampling  time  required  to  arrive  at  a  satisfactory  and  consistent 
estimation  of  sediment  particle  size  distribution  and  mass  concentration. 

This  study  represents  the  first  stage  of  characterizing  die  backscatter  amplitude 
performance  of  the  CASP,  and  implements  a  technique  to  estimate  mean  particle  size  and 
sediment  mass  concentration.  While  the  wire  target  backscatter  sensitivity  calibrations  were 
not  fully  successful,  die  sample  volume  characteristics  woe  successfully  compared  with 
theoretical  values.  Baseline  performance  data  have  been  measured  for  a  typical  sediment 
population  over  a  wide  range  of  concentrations,  enabling  the  system  sensitivities  to  be 
estimated. 


76 


APPENDIX  A 
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APPENDIX  B 


LIST  OF  SYMBOLS 


a  equivalent  particle  radius 

ag  geometric  mean  radius 

a«  transducer  radius 

c  sound  speed  in  fluid 

/«,  backscatter  form  factor 

kc  compression  wave  number 

n(a)  particle  size  spectral  density 

pa  backscattered  pressure  amplitude 

p»  on  axis  pressure  amplitude  at  reference  distance 

r  perpendicular  distance  from  transducer  to  detected  volume 

r*  reference  distance 

D  transducer  directivity 

J,  cylindrical  Bessel  function  of  order  1 

Kc  digital  counts  of  backscatter  amplitude 

M  mass  concentration 

S  overall  system  sensitivity 
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sensitivity  calculated  from  wire  targets 
voltage  output  of  receiver 
attenuation  coefficient  in  water 
attenuation  coefficient  of  scatterers 
spreading  angle  of  transducer 
half  power  angle  of  the  transducer 
Neumann  factor 
phase  shift  of  nth  wave 
acoustic  wavelength 
density  of  scatteicr 

variance  of  lognormal  distribution  about  geometric 

duration  of  transmitted  pulse 

combined  attenuation  by  fluid  and  scatterar 
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